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A stochastic difTerential equation for the plasma density dynamics is derived, consistent with 
the experimentally measured distribution and the theoretical quadratic nonlinearity. The plasma 
density is driven by a multiplicative Wiener process and evolves on the turbulence correlation time 
scale, while the linear growth is quadratically damped by the fluctuation level. The sensitivity of 
intermittency to the nonlinear dynamics is investigated by analyzing the Langevin representation of 
two intermittent distributions, showing the agreement between the quadratic nonlinearity and the 
gamma distribution. 
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The scrape-off layer (SOL) plasma of magnetic con- 
finement fusion devices exhibits a universal intermittent 
behavior characterized by a strongly non-Gaussian statis- 
tics and spatiotemporal correlations [1]. The intermittent 
character of fluctuations is closely related to a bursty 
convective transport of over-density coherent structures 
(blobs) and edge localized mode filaments (ELMs) [2]. 
Basically the deviation from normality is caused by the 
quadratically nonlinear turbulence , coupling the plasma 
density to the electric potential, i.e. + • • • cx; V"/" [3]. 
For instance an assumed Gaussian initial condition ■0o 
evolve to a Chi-square random variable (r.v.) (-^g) and 
so on. 

The closure procedure reducing the coupled determin- 
istic turbulence equations to a nonlinear Langevin equa- 
tion is still a challenging problem [4] , although a substan- 
tial theoretical effort has been devoted to the statistical 
characterization of intermittency in a turbulent media. 
In climatology science [5] as well as in fusion plasma re- 
search [6], an attempt to the dynamics interpretation of 
intermittency and its statistical signature have been ad- 
dressed by describing the intermittent quantity -0 as a 
quadratic polynomial of a Gaussian variable i.e. 

m = 9{t)+^9\t), (1) 

where lo is the non-normality parameter which mea- 
sures the deviation from the Gaussian statistics and 
the strength of the nonlinear coupling. It is not sur- 
prising that the process Eq. (1) satisfies the universal 
parabolic relation between the Kurtosis and the Skew- 
ncss, i.e. K = aS^ + b, observed in several turbulent 
media [4] , because the variable -0 -I- l/4a; is distributed as 

u! (^g + l/2w) , which is the non-central Chi-square dis- 
tribution by construction, and is a particular case of the 
gamma distribution measured in the edge of fusion de- 
vices and satisfying K = 1.55^ -I- 3 [7, 8]. Unfortunately 
the underlying physical mechanism of Eq. (1) is still not 
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clear and the Gaussianity assumption of the dynamical 
variable g is too strong. 

Based on a physical intuitions, several other existing 
stochastic models arc able to explain the emergence of the 
gamma statistics from a turbulent media. In their inves- 
tigation of scattered radiation from a fiuctuating back- 
ground, Jakeman et al. [9] showed that the gamma prob- 
ability distribution function (PDF) appears as a limit dis- 
tribution of a finite sum of independent and identically 
distributed random perturbations Xi, 

N 

^(i)=E^- (2) 

where their number N obeys a birth-death-immigration 
process with a respective rates b, d and m [9]. Without 
any condition on the perturbcrs distribution P(xi), and 
when the death rate is close to the birth rate (6 ~ d), 
then X is gamma distributed with /x = (x) scale factor 
and V = m/b shape factor. 

On the other hand, in Rcf. [10] the gamma distribution 
is derived by applying the Campbell's theorem [11] to the 
plasma density signals, assumed to be a linear superpo- 
sition of K bursts, i.e., 

K 

m = Y.^'^F{t-u) (4) 

1=1 

with exponentially distributed intensity P{x) = 
l//iexp(— a;//i), waiting time between bursts arrivals 
P{t) = cxp{—t/Tyj) and burst life time F{t) = 
exp(-t/rd). 

It is noteworthy that both stochastic processes given 
by Eq. (2) and Eq. (4) lead exactly to the same PDF 
Eq. (3), when the waiting time and the duration time 
are given by r^, = m^^ and Td = b^^ ~ d^^ . 
The equivalence between both derivations of the gamma 
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statistics has a deep physical meaning, since we could 
make correspondence between the immigration process 
and the waiting time, and between the birth process 
and the duration time. The similarities between these 
two gamma processes is extended beyond the univariate 
distribution by investigating their temporal correlation. 
The covariance function of the process Eq. (4) is calcu- 
lated using a random noise properties [11], and assuming 
an exponential burst life time distribution. 



at) 



exp{-t/Td), 



(5) 



showing that the correlations are introduced only 
through a single burst duration, since a given burst is 
independent of each others as in the Kubo- Anderson pro- 
cess [12]. The correlation structure of the shot noise pro- 
cess Eq. (4) is consistent with that of Eq. (2), which is 
also exponential with the inverse death rate as a correla- 
tion time C{t) cx exp(— &t) [13]. 

The two models provide a simple physical picture of 
edge plasma intermittency, where blobs are born close to 
the last closed flux surface (LCFS), immigrate across the 
SOL convected by the cross-field velocity, before disap- 
pearing through dissipation and parallel transport. Al- 
though the model Eq. (4) is consistent with the experi- 
mental measurements, it obeys a linear stochastic differ- 
ential equation [14], what is surprising because the in- 
termittent statistics in turbulent plasma is often associ- 
ated to a nonlinear turbulence models like the Hasegawa- 
Wakatani system equations [3]. Indeed the stochastic 
models Eq. (2) and Eq. (4) could explain the observed 
gamma statistics from the response function of the in- 
strumental devices point of view e.g., Langmuir probe, 
but their dynamical content is still far from the theoret- 
ical predictions. 

In this Letter, our purpose is to provide an insight 
of the bridge between the nonlinear dynamics and the 
intermittent statistics in the edge plasma of magnetic fu- 
sion devices. Here the nonlinear dynamics is investigated 
starting from the experimentally measured distribution. 
Namely, we assume the gamma statistics together with a 
quadratic nonlinearity as a theoretical constraint to de- 
rive a two parameters plasma turbulence equation for a 
given correlation time and fluctuation amplitude 



R= (^((V^- W)'))'^'/W = v^, 



(6) 



where < ... > references to the time average. Indeed the 
gamma distribution Eq. (3) obeys the Pearson equation, 
given in its general form by d^,P{il}) = F{ip)P{ip)/ H{iJ;). 
The associated stochastic process is specified by the fol- 
lowing Fokker-Planck equation for P{ip,t) [15], 
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FIG. 1. a quiet time series of the fluctuating plasma density 
normalized to its local mean value, with R = 30% and Tc — 3 
/is, as for a typical plasma conditions near the LCFS. 



gamma distribution Eq. (3), 



dPW (1 - l/iy) (iP) - ^ 



dip 



{iP) ip/i 



(8) 



where (■0) = z^/i being the average value of the plasma 
density field. Using Eq. (7), we derive the corresponding 
Fokker-Planck equation 



dP 
dt 



d_ 

dip 
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(9) 

where we have introduced a characteristic time scale 
r which is specified farther. Then ip follows the Cox- 
IngersoU-Ross process [16], 

dt 



dtPit) = {{iP) - m) - + v/2 Wm/'^r dw{t), (10) 

where w is the Wiener process. The stochastic model 
Eq. (10) follows a root square nonlinearity because of the 
term oc and is broadly used in financial forecasting 
[17]. The physical mechanism leading to a such dynam- 
ics is not trivial, since the typical turbulence models are 
quadratically nonlinear [3]. Hence the irreducible rep- 
resentation of the gamma distribution Eq. (8) gives rise 
to some nonlinearity, but it is still not sufficient to cap- 
ture the quadratically nonlinear dynamics expected by 
the plasma turbulence equations. However, Eq. (8) is 
degenerate and permits one to derive a higher order non- 
linear stochastic differential equation with the gamma 
PDF as a marginal. By multiplying the denominator 
and the numerator in the right hand side of Eq. (8) by tp, 
we obtain the new Pearson representation of the gamma 
process, 



dPiip) (1 - l/i^)V'" V'V(^> 



d-p 



P{^), 



(11) 



One get the following Pearson representation of the 



then the corresponding stochastic differential equation 
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follows from Eq. (7), 



(12) 



with 7 = (l + _R^)/r. In order to clarify the role of the 
time scale r, the dynamics of the correlation function p{t) 
is required. It is straightforwardly derived from Eq. (12), 



dpit) 7 
ipit) = ^ 



dt 



4'\mo)), (13) 



the term between brackets is a two-time correlation and 
can be approximated by, 



^2(t)^(0)) =r?p(t) 



(14) 



where the constant coefficients are fixed using the two 
particular cases Tc = 0, p{t) = and Tc = oo, p{t) = 1. 
We get K = {t]]) (ij^) and 77 = (V^^^ _ (^^'^ (^^2y r^j^g^^ 

Eq. (13) reduces to. 



dpjt) 
dt 



+ ip{t) = 0, 



(15) 



the correlation function is exponential pit) = exp(— 7t), 
with Tc = 1/7 as a correlation time. Then r = (1 + R^)tc 
is nothing but the correlation time increased by a fac- 
tor 1 + R^. It is worth noticing that the approximation 
Eq. (14) becomes exact in our case where -0 is gamma 
distributed. This point is not developed here and can be 
demonstrated by writing ip as a, square of Gaussian r.v. 
and using the Wick factorization theorem. 

Equation (12) is the main result of this Letter. It is 
quadratically nonlinear and evolves on a turbulence cor- 
relation time scale. The process Eq. (12) has a consistent 
physical meaning, since the linear growth is driven by the 
fluctuations amplitude. As a response to this growth, 
a nonlinear damping takes place and is proportional to 
the fluctuation level ip{t)/ The dynamics is driven 
by a multiplicative white noise with variance R^. This 
stochastic drift disappears and the dynamics evolves de- 
terministically when the fluctuations amplitude vanishes 
(i? — 0). Unlike Eq. (1) which is based on the Gaussian 
assumption of the dynamical variable, and the process 
Eq. (4) obeying a linear stochastic differential equation, 
our stochastic model Eq. (12), shows how the gamma 
distribution and the associated universal K-S scaling rise 
self-consistently from the quadratically nonlinear dynam- 
ics of the intermittent variable. Equation (12) behaves 
like the stochastic generalization of the L-H transient 
paradigm equation [18], with the difference that the lat- 
ter is quadratic for the energy E = (('(/' — (ip))/ (ip))'^, 
leading to the Nakagami distribution for the plasma den- 
sity. The stationary solution of Eq. (12) is the gamma 
distribution with stationary but not necessarily homoge- 
neous parameters. Hence it can be used to analyze the 
nonlinear time series of the plasma density at different ra- 
dial SOL positions x, using the local average field (1/') (x). 
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FIG. 2. a weakly intermittent time series of the fluctuating 
plasma density normalized to its local mean value, with R = 
50% and Tc — 16 ps, as for a typical plasma conditions of the 
near SOL. 
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FIG. 3. a strongly intermittent time series of the fluctuat- 
ing plasma density normalized to its local mean value, with 
R = 90% and Tc = 55 ps, as for a typical far SOL plasma 
conditions. 



correlation time Tc{x) as well as R{x). In order to simu- 
late a plasma density time series in different conditions, 
we have used the Ito's representation of Eq. (12) for the 
fiuctuation level ^ = ip{t)/ 



dt- 



e 



Tc(l+i?2) 



dt + R. 



Tc(l+i?2) 



^dw{t). (16) 



In Fig. 1 is plotted a quiet time series of ^(t) with time 
step dt ^ 1 ps, the fluctuations amplitude is of i? = 30% 
and the correlation time of = 3 /iS, as in the typi- 
cal edge plasma conditions close to the LCFS [7]. Fig. 2 
gives another time series representative of an intermedi- 
ate situation corresponding to the near SOL with mod- 
erate fluctuations amplitude R = 50% and Tc = 16 ps, 
tree bursts exceeding the average value by a factor of 3 
rise in the range of 1 ms. A strongly intermittent time 
series is plotted in Fig. 3 with R = 90% and r^, — 55 ps, 
tree bursts exceeding 4 (ip) and a super burst with the 



4 



amplitude of 10 (■0) are observed in the time range of 1 
ms, as is typically the case in the far SOL [6]. 

In order to investigate the role of the nonlinear dynam- 
ics in the statistical characterization of intermittent fluc- 
tuations, we compare the gamma process to the existing 
stochastic models. The beta distribution has been used 
to fit the plasma density data in the edge of TORPEX, 
according to its capability to have negative skcwness and 
its finite support [8]. Let consider here the standardized 
beta distribution, 

P(V') oc (1 - i^f'^ , < V < 1, a, ^ > 0, (17) 

where a ^ l/R^ - (-0) (1 + l/R^) and /? = a(l/ - 
1). The beta distribution obeys the Pearson equation, 
then using Eq. (7), we derive the following stochastic 
differential equation, 

7^- 4' + - ^' + ^\/2 -w{t), 18 

at T T V T 

showing that the quadratically nonlinear dynamics of 
the beta process is systematically accompanied with a 
root square nonlincarity. Thus unlike the beta process, 
the gamma process and its stochastic representation 
Eq. (12) is in good agreement with the theoretical 
quadratic nonlinearity prediction and the experimental 
measurements. Indeed it has been pointed out that 
it is sufficient to fit the plasma density signals using 
a limiting beta distribution which coincides with the 
gamma PDF for positive skewness and using a shifted 
gamma r.v., 2 (-0) — 0, for negative skewness [8]. The 
comparison between the gamma and the beta processes 
shows the importance of the PDF functional form 
used to fit experimental data. With identical two first 



moments, different distributions (gamma, beta and log- 
normal) could provide a reasonable fit of experimental 
measurements, although the nonlinearity degrees of 
their dynamics is different. Therefore, and in order to 
improve the connection between the theory of nonlinear 
dynamics and the experimental time series, the analysis 
of plasma's fluctuations using strong criteria [19-21] 
are suitable to clarify how long the statistics follows a 
given distribution. This procedure is preferable to the 
data fitting and a statistical signature based on the K-S 
scaling, since no unicity exists between this scaling and 
the underlying distributions. 

As a summary, we have derived a quadratically nonlin- 
ear stochastic differential equation for the intermittent 
plasma density in the SOL of fusion devices. The plasma 
density dynamics evolves on a turbulence correlation 
time scale (tc = 1 — 100 /is), and is characterized by 
the local fluctuations amplitude i? = 10 — 90%. This 
equation is consistent with the experimental measure- 
ments and a theoretical prediction, since it behaves a 
gamma marginal and a quadratic nonlinearity. Its simple 
representation is suitable for edge plasma modeling, 
since a flnitc plasma density correlation time, together 
with high fluctuations amplitude affect considerably the 
plasma-wall interaction and the associated transport 
of sputtered impurities and released molecules [22, 23]. 
The sensitivity of the nonlinearity degrees to the inter- 
mittent distribution has been investigated by comparing 
two potentially used distributions. Our result shows 
that the gamma distribution agrees with the quadratic 
nonlinearity, in comparison with the beta distribution 
leading to a root square nonlinearity. 
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